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Abstract 

^1^1 We examine in detail radiative corrections to the lightest scalar Higgs boson 

mass due to the top quark and scalar quark loops in the next-to-minimal super- 
symmetric standard model (NMSSM). We take into account the nondegenerate 
state for the top scalar quark masses. In our analysis, the mass matrix of the 
^ ' top scalar quark contains the gauge terms. Therefore our formula for the scalar 



. Higgs boson mass matrix at the 1-loop level includes the contribution of the 

gauge sector as well as the effect of the top scalar quark mass splitting. Thus 
we calculate the upper bound on the lightest scalar Higgs boson mass using our 
formula. We find that the absolute upper bound on the 1-loop corrected mass of 
the lightest scalar Higgs boson is about 156 GeV. 
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I. INTRODUCTION 



One of the main motivations for a supersymmetric extension of the standard model 
(SM) is the fact that a calculation of the 1-loop corrections to the Higgs boson mass 
yields a quadratic divergence arising from the SM particle loops. In a supersymmetric 
theory, all particles in the SM are accompanied by their superpartners. Therefore, all 
quadratic divergences are eliminated by the cancellation between ordinary particle and 
its superpartner loops. 

As a matter of fact, supersymmetry (SUSY) requires the existance of supermulti- 
plets made up of fermions and bosons with equal masses. If SUSY is to be relevant for 
the physical world, it must be broken, either softly or spontaneously. That is, SUSY 
must broken because in experiment one does not observe degenerate Bose-Fermi mul- 
tiplets. Thus the Higgs boson mass at the tree level receives quadratic corrections that 
are limited by an incomplete cancelation of the SM particle and its superpartner loops. 

In the SM, only one Higgs doublet is required to give masses to the quarks and 
leptons. SUSY requires at least two Higgs doublets. One is needed to give masses to 
the up-type quarks and the other to the down-type quarks and leptons. This is the 
minimal supersymmetric extension of the SM. The minimal supersymmetric standard 
model (MSSM) is the most widely studied supersymmetric extension of the SM. The 
physical mass spectra of the Higgs bosons consist of two neutral scalars, one neutral 
pseudoscalar, and a pair of charged scalars. According to the tree level potential, the 
mass of the hghtest scalar Higgs boson have to be lighter than the Z boson mass, and 
the the charged Higgs boson mass must be heavier than the W boson mass. 

Recently radiative corrections to the Higgs boson masses have been calculated by 
many authors. The tree level constraints on the charged Higgs boson mass can be 
violated when the top quark mass is heavy and the pseudoscalar Higgs boson mass is 
light. This means that radiative corrections to the charged Higgs boson mass can be 
negative contributions. However, these corrections to the tree level result are numer- 
ically very small [1,2]. On the other hand radiative corrections to the lightest scalar 
Higgs boson mass give a significant contribution and the mass of the lightest scalar 
Higgs boson can be substantially heavier than the Z boson mass [3-5] . 

If one considers about more general extensions of the Higgs sector [6] , the tree level 
results in the MSSM can be changed. The simplest extension is the next-to-minimal 
supersymmetric standard model (NMSSM). The Higgs sector in the NMSSM consists 
of two Higgs doublits Hi and H2 and a Higgs singlet chiral superfield A^. In this 
model, the superpotential contains a new coupling among two Higgs doublits and a 
Higgs singlet. Thus the parameter in the superpontial of the MSSM can be generated 
dynamically by the vacuum expectation value ( VEV) of the singlet Higgs field [7] . The 
NMSSM has ten real degrees of freedom. After the Higgs mechanism takes place, three 
degrees of freedom correspond to neutral ( G ) and charged ( G^, ) unphysical 
Goldstone bosons, the other seven correspond to neutral ( 5*1, 5*2, 6*3, Pi, P2 ) and 
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charged ( C"*", C~ ) physical particles. Here 5*1, 5*2, and 5*3 are scalar Higgs bosons 
while Pi and P2 are pseudoscalar Higgs bosons. As a result, the effective number of 
parameters describing the Higgs sector at the tree level is six. 

As we have already mentioned above, the tree level result in the NMSSM is a little 
different from that in the MSSM. The mass of the charged Higgs boson may be lighter 
or heavier than the W boson mass. Also the upper bound on the lightest scalar Higgs 
boson mass can be heavier than the Z boson mass [8] . Therefore the upper bound on 
the lightest scalar Higgs boson mass of the NMSSM at the tree level is heavier than 
that of the MSSM. Furthermore radiative corrections to the lightest scalar Higgs boson 
mass lead to a important contributions. It is well known that the largest contribution 
to the hghtest scalar Higgs boson mass comes from the top quark and scalar quark 
loops [9]. These effects are quite substantial when the top scalar quark mass is much 
heavier than the top quark mass. This implies that at least one top scalar quark 
mass is heavier than top quark mass. If the top quark and scalar quark masses are 
identical, the contribution of radiative corrections to the lightest scalar Higgs boson 
mass vanishes. 

Assuming the degeneracy of the left- and right-handed top scalar quarks, radiative 
corrections to the lightest scalar Higgs boson mass are studied in ref. [10,11]. Especially, 
the anthers of ref. [11] performed a low energy renormalization group analysis of the 
Higgs sector of the model, and they arrived at the conclusion that the upper bound on 
the lightest scalar Higgs boson mass is 123 GcV for irit = 180 GeV if the masses of the 
top scalar quark were assumed to be degenerate at Msusy = 1 TeV. 

The main aim of this paper is to explore all the parameter space and provide 
an accurate numerical evaluation for the upper bound on the lightest scalar Higgs 
boson mass in the NMSSM. We examine in detail the radiative corrections to the 
hghtest scalar Higgs boson mass using the 1-loop effective potential. In general, the 
contribution of radiative corrections due to the bottom quark and scalar quark loops 
to the tree level mass of the lightest scalar Higgs boson is very small numerically. 
Therefore we use the effective potential which includes the contribution of the top 
quark and scalar quark loops assuming the nondegenerate of the top scalar quark 
masses. Thus the mass splitting among the top scalar quark give a significant effect to 
the upper bound of the lightest scalar Higgs boson mass. 

In fact, there have been many studies about the upper bound on the lightest scalar 
Higgs boson mass including radiative corrections in the NMSSM. Unhke their calcula- 
tions we especially use the top scalar quark mass matrix containing the gauge terms 
in the 1-loop effective potential. Therefore our formalism for the scalar Higgs boson 
mass matrix at the 1-loop level contains the contributions of the gauge sector. After 
including radiative corrections, the upper bound on the lightest scalar Higgs boson 
mass becomes dependent on mg, my. At, and rut as free parameters as well as the pa- 
rameters at the tree level. Of course our formalism for the mass matrix is given by the 
complicated expression. Nevertherless we believe that the upper bound on the lightest 
scalar Higgs boson mass obtained by our formalism can provide a reliable result. 
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II. HIGGS BOSON MASSES 



The scalar Higgs potential of the NMSSM comes from the auxiliary F- and D-fields 
and the soft SUSY breaking terms. The Higgs boson masses in the model at the tree 
level are given by the relevant potential 

V^^Vp + VD + Vsoft (1) 

which is expressed in terms of two Higgs doublet fields Hi and H2 and a complex 
singlet field N as 

-lxk*HiH2N*^ + H.C.) , 

Vn = 4(Hl^H2 + HlaHir + ^-l(\H2\'-mr , (2) 

o o 

- i^kAkN' + H.c.) . 

In the above expressions, gi and g2 are the U(l) and SU(2) gauge coupling constants, 
respectively, a — {a^,a'^,a^) are the Pauli matrices, Ak and Ax are the trilinear soft 
SUSY breaking parameters, and H1H2 — i??-f^2 ~ H~H'^. The parameters A and k 
come from the trilinear couphngs in the part of the superpotential of the NMSSM 

XNH1H2 - -N^ . 
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We assume that only neutral components of the three Higgs fields Hi, H2, and 
N acquire vacuum expectation values Vi, V2, and x, respectively. Without loss of 
generality, we also assume that the vacuum expectation values are real and positive. 
The conditions that the Higgs potential becomes minimum at Hi — vi, H2 — V2, and 
N — X yield three constraints, which can eliminate the soft SUSY breaking parameters 
'iTT'Hi, fnH2-i ^iid mjv from the potential Vq. 

Then the mass spectra of the Higgs bosons can be expressed in terms of six pa- 
rameters, A, k, Ax, Ak, X, and tan/3 = V2/V1. Note that m| = {gi + g^v^ {2 for 

V = \J I'l + t'2 = 175 GeV. Collecting the real parts of the neutral Higgs fields from Vq, 
one can obtain the mass spectrum of the scalar Higgs bosons at the tree level. It is 
given by a symmetric mass matrix, whose elements are as follows: 

M^i = (mz cos (3f + Xx tan (3{Ax + kx) , 

M^2 = {mz sin + Xx cot l3{Ax + kx) , 

M°3 = {2kxf - kxAk + i^v'^^x sin 2(3 , (3) 

= (A^^;^ - ^m|) sin 2(3 - Xx{Ax + kx) , 
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M2°3 = XV sin (3 — \v cos I3{A\ + 2A;x) , 
cos j3 — \v sin /3(74;^ + 2kx) . 

When the above mass matrix is diagonahzed, we would obtain squared masses for the 
three scalar Higgs bosons. The smallest eigenvalue yields the mass of the lightest scalar 
Higgs boson, ms^- 

The tree-level masses obtained from Vq are subject to radiative corrections. These 
corrections are contributed in supersymmetric models by the quark and scalar quark 
loops. According to the Coleman- Weinberg mechanism [12], the 1-loop corrections are 
given by the effective potential containing the mass squared matrix of the scalar 
quark, 

where the supertrace is defined as 

Strf{M') = E(-l)''"(2^. + . (5) 

i 

The squared masses of the quark and scalar quark are denoted by ml and their spins 
are Jj = for the quark and Jj = 1/2 for the scalar quark, respectively. The arbitrary 
scale A is taken to be Msusy = 1 TcV. Thus the expression for the 1-loop effective 
potential which corresponds to taking into account the radiative corrections due to the 
top quark and scalar quark loops is given as 

- 32^<i^°SA^-2j + 32^<i^°SA^-2j 

7:rnl log —7 , 6 

167r2 « \^ ^ A2 2; ' ^ ' 

where ruq^ and m^-j are the masses of the top scalar quark, ruq^ < rriq^ and rrit is the 
top quark mass. 

In terms of Vi, the mass matrix of the scalar Higgs bosons including radiative 
corrections at the 1-loop level is given by 

, _ / d'V, _ \_d_V, \ 

where 0j are the conventional notations for the real and imaginary parts of the Higgs 
fields. The first derivative term of the effective potential is minimization condition. 

In order to calculate the radiative corrections, one has to know the squared masses 
of the scalar quarks in V\. They are given as the eigenvalues of a 4 x 4 Hcrmitian matrix 
for the third generation. The scalar quark mass matrix may be broken up into four 
2x2 submatrices, because only two block- diagonal submatrices depend on the neutral 



(7) 

H'^=v^.H'^=vo.N=x 
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Higgs sector and thus arc nonzero. The two eigenvalues of the upper-left block-diagonal 
submatrix are the squared masses of the left- and right-handed top scalar quarks, and 
those of the lower-right block- diagonal submatrix the squared masses of the left- and 
right-handed bottom scalar quarks. Although the contribution of the bottom quark 
and scalar quark loops for large tan (5 is not negligible, generally the leading radiative 
corrections to the scalar Higgs bosons are the contribution of the top quark and scalar 
quark loops. 

Taking only the top quark and scalar quark into account, the matrix elements of 
Ai^ depending on the neutral Higgs sector are: 
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Mn = m'Q + hf\Hl 



— (l^iT 



(8) 



where ht is the Yukawa coupling of top quark, in the part of the superpotential of the 
NMSSM 

htQf^H^ , 

and TTiQ and rriT are the soft SUSY breaking scalar quark masses and is the charge 
conjugate of the right-handed top quark. We assume that the soft SUSY breaking 
scalar quark masses are real and positive. 

The eigenvalues of the matrix are expressed by the top scalar quark masses 



m 



tl,2 



1 1 

+ 2("^0 + ^t) + ^ml cos 2/3 



± 



1 2 5 \ ^ 

-(m^-m^) + (-m^- — m|)cos2/3j + ml{At + Xx cot pf 



(9) 



htV sin f3 and ml 



fi'K^i + '^2)/2- From the formula we find that the 



where rrit — lotvauLfj ain_i llb^, 
masses of the top scalar quark are symmetric under interchange mg 



rriT- 



Substituting these expressions into Vi, one can now calculate the radiative correc- 
tions. After lengthy calculation, the elements of the mass matrix for the scalar Higgs 
bosons including the 1-loop corrections are obtained as 



Ml, = 



87r2 
+ 



2 I 2 2^ , 3m|cos^/3 ^ rn}m} 



2mfAtXx 
sin 2(3 



A4 



+ 



m| cos P 



'3 
Ai 



-m^ - -ruz) cos /3 



2 2 

m:- — mf- 



log 



m 



ti 



2 

mf- 

12 
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1 ^ 

in^v^ sin^ p A2 



3mf 



m^AtXx cot f3 4 



TT — m| sin /3 



4 2 2\2-2^ 




^ — 7T - 7:"^7 sm B log — 7-r 



327r2^;2 



3 

— m^A^ cot^ /3( A + Ax cot I3fg{m\ , m| ) 



+ 



m 



' AtX cot Pf{ 



87r2 
+ 



+ 



327r2v2 
3 

32^ 



-m^- -m|)2sin2/3- * 

D 



4m? 



m| cos /3A2 + ( — - m|) sin /3A 

sin p 




87r2 



A2m^Acot/5(74t + Ax cot P)g{m}^,m}^ 



+ 



3m^X cot P / 4m: 
32^^2 (^—^ 

_3m|AAcot^ 2 2^ 
167r2t;sin/3 ^""fi ' ""t^ i 




87r2 
+ 



Aim(Acot/9(At + Ax cot P)g{m\,m\) 

SmimfX cos /3 cot /? / At + Ax cot /9\ 
Tz::-^ n n — loe 



m 



'ti 



32772^; 



2 2 



m 



t(A + 2Ax cot P)f(mf , mj- ) 

167r2't;sin/3^ * ^ 



'f2 



\ A At + Ax cot/?) 
sm /? 

cos/3 [ 2 2n ^ 

+ Axcot/3) 



T;sin P 

sin /? 



A 2 



(m^ - m^) + (; 



■m 



w 



^m|)cos 2p 



— g^l) cos 2P 



A 2 



The two functions / and g are defined as 



f{ml,ml) 



g{ml,ml) 



(mi - ml) 
1 



{ml — m|)^ 



1 




(12) 



Note that these calculations have been done by several authors, either with some 
approximations or without [13-16]. Nevertheless, their formulas for the mass matrix of 
the scalar Higgs boson have not contain the contributions of the gauge sector at the 
1-loop level. Unlike their calculation we have derived a formula for the mass matrix 
of the scalar Higgs boson including the contribution of the gauge sector. Furthermore, 
the number of the parameter for the mass matrix at the 1-loop level is the same that 
of their formula. When radiative corrections to the tree level mass of the scalar Higgs 
boson are included, the values of the top quark Yukawa coupling are determined by the 
equation, mt = htv sin f3. As a consequence, one could have taken the top quark mass 
instead of the top quark Yukawa coupling as a parameter. Thus the mass matrix at the 
1-loop level contains mg, mx, A^, and mt as free parameters as well as the parameters 
at the tree level. 



Now let us consider the mass splitting among the top scalar quark. After consider- 
ing all conditions and assumption in the NMSSM, we initiate a numerical study for the 
mass spliting among the top scalar quark. In the renormalization group analysis of the 
NMSSM, the maximum values of the dimensionless parameters A and k are 0.87 and 
0.63, respectively. Prom the experimental result of the CDF and DO collaborations, the 
range of top quark mass is given as 176±8±10 and 1991^? ±22 GeV, respectively [17]. 
Therefore, we take the top quark mass to be 175 GeV. The range for tan/9 is defined 
in terms of many ways. It is demonstrated that the renormalization group analysis of 
radiative symmetry breaking always leads to tan (3 values that are larger than 1 [8] . In 
our argumentation, we consider the ratio of VEV's tan/3 = V2/V1 to be in the range 



Thus the maximum value of the free parameter tan /3 is about 39 for mt — 175 and m^ 
= 4.5 GeV. From the values of tan/3 and mt, the value of the coupling A is given by 
the analysis of the renormalization group equation. 

Here we assume that mf^ is heavier than the mass of the top quark. For the 
convenience we then define the ratio of the top scalar quark masses by R, 



III. NUMERICAL ANALYSIS 



2 < tan/3 < — . 

rub 



(13) 



< R 



< 1 . 



(14) 



8 



Thus R is the inverse of the mass sphtting among the top scalar quark. This imphes that 
the mass sphtting among the top scalar quark increases with decreasing R values, vice 
versa. The curves for the maximum and minimum of R are denotted by Rmaximum and 
Rminimum in Fig. 1. As Can be seen in the figure, these curves decreases monotonically 
as the value of At increases to 3000 GeV. That is, the mass sphtting among the top 
scalar quark always increases with increasing At value. 

Next we now wish to evaluate the numerical result for the upper bound on the 
lightest scalar Higgs boson mass including radiative corrections in the NMSSM. The 
radiatively corrected squared masses of the scalar Higgs bosons are obtained as the 
eigenvalues of M° + M^. It is not easy to obtain an analytic expression for the smallest 
eigenvalue, in terms of the various parameters that are present in + M^. 

The hghtest scalar Higgs boson mass at the 1-loop level is obtained by the numerical 
diagonalizion of the mass matrix of the scalar Higgs boson. 

Fig. 2 shows the upper bounds on the lightest scalar Higgs boson mass including 
the radiative corrections due to the top quark and scalar quark loops, as a function 
of X, for mt = 175 GeV and tan/? = 2, < A < 0.87, < k < 0.63, and 100 GeV 
< Ax{Ak,mQ,mT) < 1000 GeV. The bounds on ms, for At = 1000, 2000, and 3000 
GeV have maximum values as the value of x approach 1000 GeV. The bound for At 
— 2000 GeV is greater than those for the other values of At when the value of x is 
given by 1000 GeV. We can infer that the bound on ms^ for tan P — 2 and x — 1000 
GeV is not increase thought the value of At is greater than 3000 GeV. Therefore the 
upper bound on ms^ for tan (3 = 2 and x = 1000 GeV is determined in terms of a value 
between 1000 and 3000 GeV. Thus we come to the conclusion that the upper bound 
on TTis^ for tan/3 = 2 and x = 1000 GeV can be greater than 143 GeV in an energy 
between At = 1000 and 3000 GeV. In Fig. 3 we plot the upper bound on ms^ for tan /3 
= 6 when the other parameters are given by the same parameter space as those of Fig. 
2. We find from Fig. 3 that the upper bound on m^^ for x = 300 GeV and At = 2000 
GeV is about 151 GeV. The value of ms^ = 151 GeV is greater than that for tan/? = 2 
in Fig. 2. Therefore, let us calculate the upper bound on ms^ for tan /? = 10 through 
the same parameter space as those in Fig. 3. Fig. 4 shows that the upper bound on 
for X = 160 GeV and At = 2000 GeV is about 148 GeV. But the value of ms, = 
148 GeV is smaller than that for tan/3 = 6 in Fig. 3. From Fig. 2, Fig. 3, and Fig. 4, 
we surmise that the upper bound on ms^ for tan /3 > 10 is smaller than that on 7715^ 
for tan /3 = 6. 

In Fig. 5 we plot the upper bound on ms^ at the 1-loop level, as a function of At, 
for mt = 175 GeV, < A < 0.87, < A; < 0.63, and 100 GeV < AxiAk,mQ,mT, x) < 
1000 GeV. The upper bound on ms^ for tan /3 = 6 is not increase though the value of 
At is greater than 3000 GeV. Again we conject that have a maximum value in a 
value between tan (3 = 2 and 10. Fig. 6 shows that the upper bounds on ms^ at the 
tree level and 1-loop level, as a function of tan/3, for rrit = 175 GeV, < A < 0.87, 
< A; < 0.63, 100 GeV < AxiAk,mQ,mT, x) < 1000 GeV, and GeV < A < 3000 
GeV. The absolute upper bound on ms^ at the tree level and 1-loop level are 101 and 
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156 GeV for tan (3 f» 2.7, respectively. 



IV. RESULTS AND CONCLUSIONS 

In order to devise effective search strategies for the detection of Higgs particle, 
the study for the lightest scalar Higgs boson mass is very important. Particularly, 
supersymmetric models impose strong constraints on the lightest Higgs boson mass. 
In the NMSSM the upper bound on at the tree level is about 130 GeV for all the 
parameter space. By the renormalization group analysis the improved upper bound on 

is about 101 GeV for rrit ~ 175 GeV. 

We know from the previous papers that the top quark and scalar quark loops have 
a positive effect on the lightest Higgs boson mass when radiative corrections to the 
Higgs sector of the tree level are considered. But authors of the papers calculated the 
lightest scalar Higgs boson mass with some approximations. Here we examine in detail 
radiative corrections to the lightest Higgs boson mass without some approximations. 

The input parameters at the 1-loop level as well as at the tree level are varied 
independently. We find from Fig.l that At play an important role in the mass splitting 
among the top scalar quark. Here the value of At is allowed to vary from GeV to 
3000 GcV. The free parameters, Ax, A^, mg, rriT, and x are limited as masses below 
1000 GeV. In this case, the upper limits on the top scalar quark masses for = 175 
GeV are approximately given by ~ 1000 GeV. We are investigated the dependence of 
the lightest scalar Higgs boson mass on all input parameters. We calculated radiative 
corrections to the lightest Higgs boson mass using the 1-loop effective potential in the 
NMSSM. Here we consider the contribution of the top quark and scalar quark loops 
for the effective potential. 

We take the mass matrix of the top scalar quark including the gauge term. Fur- 
thermore the mass matrix of the top scalar quark arc given as the nongeneratc state. 
Additionally we assume that the lower limits on the top scalar quark masses are greater 
than the top quark mass. Thus radiative correctons to the lightest Higgs boson mass 
contain the effect of the mass spliting among the top scalar quark and the contribu- 
tion of the gauge sector. An analytic formula for the scalar Higgs boson mass matrix 
including these effects are derived. The upper bound on m^^ containing these effects 
is calculated by our formula. 

We find from Fig. 5 that the upper bound on m^-^ have a maximum value as At 
approach a particular value. But the value of At can not exceed 3000 GeV when the 
upper bound on ms^ is maximised. In other to find the absolute upper bound as a 
function of tan/3 including radiative corrections, we maximized over the parameter 
space using numerical analysis. The result of our calculation are presented in Fig. 6. 
We also find from Fig. 6 that there exists the absolute upper bound on 777.5^ of about 
156 GeV. 
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Figure Captions 



Fig. 1. The mass ratio of the top scalar quarks, as a function of At, for rrit = 175 
GeV, < A < 0.87, 2 < tan/3 < 10, and 100 GeV < mQ(mr,x) < 1000 GeV. The 
upper and lower curves are the maximum and minimum of R, respectively. 

Fig. 2. The upper bound on the hghtest Higgs boson mass at the 1-loop level, as a 
function of x, for rut = 175 GeV, tan/3 = 2, < A < 0.87, < k < 0.63, and 100 
GeV < Ax{Ak,mQ,mT) < 1000 GeV. These curves correspond to A = 1000, 2000, 
and 3000 GeV, respectively. 

Fig. 3. The same as Fig. 2, except for tan/3 — 6. 
Fig. 4. The same as Fig. 2, except for tan/3 = 10. 

Fig. 5. The upper bound on the lightest Higgs boson mass at the 1-loop level, as 
a function of A, for mt = 175 GeV, < A < 0.87, < A; < 0.63, and 100 GcV 
< Ax{Ak,mQ,mT,x) < 1000 GeV. These curves correspond to tan/3 — 2, 6, and 10, 
respectively. 

Fig. 6. The absolute upper bound on the lightest Higgs boson mass, as a function of 
tan /3, for mt = 175 GeV, < A < 0.87, < A; < 0.63, 100 GeV < ^(Afc, tuq, mr, x) < 
1000 GeV, and GeV < At < 3000 GeV. These curves correspond to the tree level 
and 1-loop level, respectively. 
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